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We report the first experimental observation of tlie two-node tliickness dependence of tlie critical 
current in Josephson junctions with a ferromagnetic interlayer. Vanishings of the critical current 
correspond to transitions into vr-state and back into conventional 0-state. The experimental data 
allow to extract the superconducting order parameter oscillation period and the pair decay length 
in the ferromagnet. We develope a theoretical approach based on Usadel equations, which takes 
into account the spin-flip scattering. Results of numerical calculations are in good agreement with 
the experimental data. 



One of the exciting topics in studying the coexis- 
tence of superconductivity {S) and ferromagnetism [F) 
is proximity-induced sign-reversal superconductivity in a 
ferromagnet close to an SF-interface P, |3| ■ The super- 
conducting order parameter does not simply decay into 
the ferromagnet but also oscillates. An undoubted ev- 
idence of sign-reversal spatial oscillations of the super- 
conducting order parameter in a ferromagnet was the 
observation of the 7r-state in SFS Josephson junctions 

0. 0.0.0. 'TT-junctions' are weakly coupled super- 
conducting structures which demonstrate 7r-shift of the 
macroscopic phase difference in the ground state. The 
relation between the superconducting current Is and the 
phase difference c/? in a Josephson junction is described by 
a 27r -periodic function. In the simplest case of a tunnel 
barrier or a barrier, made of dirty normal metal, one finds 
Is — /cSint/?. The Josephson 7r-junction has an anoma- 
lous current-phase relation Is = /c sin((^-|-7r) = —Ic siniy9, 

1. e. it is characterized (nominally) by the negative criti- 
cal current 0- Spatial oscillations of the superconduct- 
ing order parameter in a ferromagnet close to an SF- 
interface was predicted in Ref. 0. A physical origin of 
the oscillations is the exchange splitting of the spin-up 
and spin-down electron subbands in a ferromagnet. It 
was discussed in Refs. H, IE 01 that in order to ob- 
serve manifestations of the transition into the 7r-state 
one should fabricate an SFS sandwich with the i^-layer 
thickness dp close to integer numbers of half-periods of 
the order parameter spatial oscillations Aea;/2. The pe- 
riod is Xex — 27r^i?2, where the oscillation (or "imag- 
inary" ) length S,F2 can be extracted from the complex 
coherence length £p in a ferromagnet: = -z^ — h It^- 
In the simplest case the imaginary length and the 
order parameter decay length S^pi are equal \f\: S^pi ~ 
^F2 = \/ fiD /Eex , where D is the diffusion coefhcient 
for electrons in a ferromagnet and E^x is the exchange 
energy responsible for sign-reversal superconductivity in 
a ferromagnet. Temperature changes of the coherence 



length related to the thermal energy contribution to pair- 
breaking processes were introduced in Ref. Q , in which 
temperature — 7r-transition was observed for the first 
time. The expressions for ^pi and are the following: 

_ j hD _ I hD irkBT 

~ V ^{TTkBTy+Ei^±TTkBT ~ SKTx 2E,x 

(1) 

The latter approximation corresponds to the case Eex S> 
fcsT, which is valid for experiments discussed below. 

Detailed experimental studies of the critical current 
thickness dependence for Nb — Cuq, 47^10,53 — Nb Joseph- 
son junctions has been started by us in Ref. 0|. A very 
large decay of the critical current and its sharp reentrant 
behavior for thicknesses close to 23 nm have been ob- 
served. An analysis of the experimental data and their 
comparison with the modern model described below have 
shown that the observed deep minimum is probably the 
reverse transition from the tt- into the 0-state at the F- 
layer thickness close to full oscillation period while the 
first node of the dependence has to be at the thickness 
of about 10 nm. Thus, the presented work is devoted 
to finding of the two-node behavior of the SFS junction 
critical current as well as to discussion of mechanisms 
of the strong order parameter decay in a ferromagnetic 
CuNi interlayer. 

In fact a nonmonotonic Ic{dp) dependence close to 
— TT-transition was observed for the first time in Ref. 0| 
and has been presented there as a number of Ic{T) curves 
for different thicknesses dp. Later Kontos et al Jj] for 
Nb ~ Pdo,gNiQ,i — Nb and then Sellier et al [3 for 
iV6— CitQ. 52 A"*o.48~Af6 junctions measured detailed reen- 
trant Ic{dp) curves for i^- interlayer thicknesses close to 
— TT-transition. In this work we have investigated the 
thickness dependence of the 5'i*"S' junction critical current 
density in a wide thickness range for sandwiches fabri- 
cated as described in Ref. 0. All junctions had their 
lateral sizes smaller than the Josephson length and uni- 
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FIG. 1: The F-layer thickness dependence of the critical cur- 
rent density for — Cuo.irNio.^s — Nb junctions at tem- 
perature 4.2 K. Open circles present experimental results, 
solid and dashed lines show model calculations discussed in 
the second part of the paper. 



form current distribution. To do this the junctions with 
F-layer thicknesses of less than 17 nm were made with 
the contact area 10 x 10 ^m^ and all the rest had the 
area 50 x 50 /zm^. Weakly- ferromagnetic Cmq. 47-^*0. 53- 
interlayers had the Curie temperature of about 60 K. In 
the thickness interval 8 — 28 nm we had about 6 orders 
of the critical current density change with vanishings at 
two dp values as it is presented in Fig. ^ One can see 
that undoubtly the curve demonstrates both direct — tt- 
transition and reverse transition from tt- to 0-state. In 
transition points the critical current Ic{dp) is equal to 
zero and then should formally change its sign. Since in 
real experiments we could measure only magnitude of 
the critical current, the dependence Iddp) between two 
sharp cusps is the negative (corresponding to the 7r-state) 
branch of the curve which is reflected into the positive 
region. Due to slight temperature dependence of the or- 
der parameter oscillation period in our weak ferromagnet 
(described by ^) we could pass through the transition 
points using samples with critical i^-layer thicknesses 11 
nm and 22 nm by means of temperature decrease. Tem- 
perature — TT- and TT — -transitions are presented in 
middle panels of Fig.|21 Upper and lower panels show the 
critical current temperature behavior for samples with 
close F-layer thicknesses. One can see that we lost a pos- 
sibility to detect temperature transitions changing the 
thickness only by 1 — 2nm. This implies that the temper- 
ature decrease from 9 K down to 1 if is accompanied by 
the decrease of 1 — 2 nm in the spatial oscillation period 
and by the decrease of about 0.3 nm in the oscillation 
length. In this temperature range the change of is 
about 0.2 nm as it has been estimated from I (.{dp) curves 
at different temperatures. At the same time simple eval- 
uations of ^Fi (obtained from the slope of the Iddp) 
envelope) and ^^^2 (estimated from the interval between 
two minima) show a large difference between these two 



lengths (1.3 nm and 3.5 nm, correspondingly) that can 
not be explained by the thermal contribution described 
by CJ. 

So, to carry out a theoretical analysis of the results 
obtained, we need to specify the nature of additional de- 
pairing processes that increase ^_f2 and decrease ^_fi. As 
the F layer is an alloy, a role of the magnetic scattering 
may be quite important |1, |^. Magnetic inhomogeneity 
is related above all to Ni-iich clusters H, 0| arising in 
Cui-xNix ferromagnet for x close to 0.5. In the region 
of these concentrations when the Curie temperature is 
small, we may expect that the inverse spin-flip scattering 
time hT~^ could be of the order of the average exchange 
fleld Ef^x or even larger. This circumstance strongly mod- 
ifies the proximity effect in the SF systems. A role of 
spin-orbit scattering should be neglected for the CuNi al- 
loy since it is substantial only in ferromagnets with large 
atomic number Z. To take into account the exchange 
fleld and the magnetic scattering in the framework of 
Usadel equations it is necessary just to substitute Mat- 
subara frequencies hy co ^ UJ+ lEex + Gh/T, 113, where G 
is the normal Green's function. Note that this procedure 
assumes a presence of the relatively strong uniaxial mag- 
netic anisotropy which prevents mixing of spin-up and 
spin-down Green functions |lll |. 

To have some idea about the influence of the magnetic 
scattering on the proximity effect we may start with the 
linearized Usadel equation for the anomalous Green's 
function in a ferromagnet 



Lu\ + iEexSgn (uj) + — j Ff — 



f 



0. 



(2) 



The exponentially decaying solution has the form 

Ff {x,uj > 0) ^ Acxp{-x{ki+ik2)) , (3) 

with 
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Here = y^hD/Kx and ^Fr,2 = 1/A;i,2. The anoma- 
lous Green's function F/ at w ~ ksTc gives us an idea 
about the spatial variation of the Cooper pair wave func- 
tion. In the limit of the vanishing magnetic scattering 
and ksTc « E^x the decaying (^fi) and the oscillat- 
ing (^F2) lengths are practically the same. However, 
if the spin-flip scattering time becomes relatively small 
EexTs/fi ^ 1, the decaying length could be substantially 
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FIG. 2: Temperature dependences of the SFS junction criti- 
cal current density at several f-layer thicknesses close to the 
critical ones. The dashed lines show calculation results based 
on Eq. 0. 



smaller than the oscillating length. This results in much 
stronger decrease of the critical current in SFS junctions 
with increase of the F layer thickness. 

We have seen it experimentally that the form of 
Icidp) dependence varies a little with temperature, so 
a good idea about this dependence may be already ob- 
tained from the temperature region near T^. Using 
k = ki + ik2 in the form 



k = 



uj\ + iEexSgn (w) H ) /HD 




(4) 



we can obtain (see Ref. fo'^ the case of good SF- 
interface transparency and dp ^ S^pi the following ex- 
pression for the critical current: 



{cOsdF/^P2 + i^Pl/^F2) sin dp /^F2) 



(5) 

where are taken in the limit of a; <C E^x, fi/Ts- 

Now we shall address the question of the exact thick- 
ness and temperature dependence of the critical current 
in SFS junctions. To deal with the complete set of 
the Usadel equations it is convenient to apply the usual 
parametrization of Green's functions : G/ = cos0(a;) 
and Ff = sin8(x). Then for > the Usadel equation 
is written as 



hcosQ 



sin 9 



= 0. 



(6) 



If the temperature variation of the exchange field is 
negligible at T < Tc the most direct way how the tem- 
perature could interfere is through the Matsubara fre- 



quencies. The presence of the magnetic scattering pro- 
vides another mechanism of the critical current temper- 
ature dependence - through the normal Green's function 
Gf — cos 8. The important range of the Matsubara fre- 
quencies variation is of the order of fc^Tc for supercon- 
ductivity. Then in the case of the relatively strong mag- 
netic scattering hT~^ >> ksT^ the second mechanism of 
the temperature dependence will be predominant. 

In the limit of relatively large F-layer thicknesses 
dp > £,Fi and rigid boundary conditions we may obtain 
an analytical solution of Eq.© and the expression for 
the critical current density reads 



jc[dp,T) = Re 

e^P 

with the function 



( 



T (n)q<iy.^{-qy) 



(7) 



(A/ {2TrkBT)Y 



n + 1/2 + ^{n + l/2f + {A/ {2TrkBT)f 



2 ' 



and 1 ~ — 



and y = dp/S^p, q =^ ^/2i + 2a + 2uj , where 

{i + u>)/{a + i + uj). 

In the hmit a ^ and ksTc « E^x Eq. Q coin- 
cides with that obtained previously in Ref. 13] . The 
theoretical fit of our experimental results which is based 
on Eq. (UI) is presented in Fig. ^ by the solid line and in 
Fig. m by dashed lines. Besides the dashed line in Fig. 
shows calculations made using Eq. One can see a 

good agreement obtained with the following parameters: 
Eex/ks ~ 850 K, h/Ts w 1.33 E^x, = 2.16 nm. The 
fitting also yields considerable value of 'dead' layers d^: 
2do « 4.3 nm, which do not take part in creating of the 
'sign-reversal' superconductivity. The dead layer may 
arise due to not full correspondence of the theoretical 
approach and the real system. On the other hand, other 
experiments |Q also demonstrate the existence of large 
enough (2-3 nm) nonmagnetic layers at S'F-interfaces. 

A final remark concerns the real transparency of SF- 
interfaces in our SFS sandwiches. In modern theories 
an interface transparency is characterized by a parame- 
ter jB = {RbS/ pp£,*), where Rb is interface resistance 
per unit area, S is the SFS junction area, pf is J^-layer 
resistivity and S,* = yKD/2/KkBTc. To estimate Rb and 
Pp we have carried out detailed measurements of SFS- 
junctions /^-characteristics. The upper inset in Fig. |2| 
shows that JU-characteristics are described by the ex- 
pression V — Ry'P — I^. The linear approximation 
presented in Fig. 13 has given Rb ~ 30 /xfi for junctions 
with the area of 10 x 10 nm? and pp « 62 p^l ■ cm. 
It allows to estimate following ferromagnet parameters: 
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of an extra mechanism of the order parameter decay in 
CuNi, mainly related to the strong spin-flip scattering 
on magnetic inhomogeneity. 
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FIG. 3: Resistance of S'.FS'-sandwiches normalized to the 
junction area of 10 x 10 /xm^ vs. the .F-layer thickness. In- 
sets show typical IV- and Fraunhofer {Ic{H)) dependences 
of SFS junctions with ideal fitting by well-known Josephson 
expressions. 

the electron mean free path I ^ 1 nm, the diffusion co- 
efficient -D « 5.2 cw? jc and the characteristic spatial 
scale ^* K, 9.4 nm. Values obtained determine the good 
enough transparency parameter 73 = 0.52 that confirms 
the validity of the approximation used. 

An additional breakthrough of the work is fabrication 
of TT-junctions with large enough critical current density. 
Solving of this problem has enabled detailed experimen- 
tal investigations of — 7r-transition peculiarities, reliable 
detections of second harmonic in the current-phase rela- 
tion and the — 7r-coexistence . High magnitude 
of the critical current also allows to use SFS 7r-junctions 
as stationary phase 7r-shifters in novel modifications of 
the digital and quantum logic |l7j . In proposed logic cir- 
cuits TT-junctions are connected together with ordinary 
tunnel junctions and should not introduce themselves any 
noticeable phase shift during dynamical switchings in the 
rest of the circuit. This is possible only if the 7r-junction 
critical current is much larger then critical currents of 
other junctions. The Nh — CuNi — Nb 7r-junctions are 
based on the standard niobium thin film technology so 
they can be incorporated directly into existing architec- 
tures of the superconducting electronics. 

Thus, both — TT and reverse tt — transitions have 
been detected in SFS {Nb — CuNi — Nb) junctions for 
the first time. The double-reversal thickness dependence 
of the critical current is a most striking evidence of the 
superconducting order parameter spatial oscillations in 
a ferromagnet close to SF-interface. We have also ob- 
served that the oscillation length in the ferromagnetic 
CuNi alloy is considerably larger than the the pair de- 
cay length. We have presented a theoretical description 
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